EXPLICIT NON-ALGEBRAIC LIMIT CYCLES FOR 
POLYNOMIAL SYSTEMS 



A. GASULL*, H. GIACOMINI, AND J. TORREGROSA* 

Abstract. We consider a system of the form x — P n (x, y) + xR m (x, y), y = 
Q n (x,y) + yR m (x,y), where P n (x,y), Q n (x,y) and R m (x,y) are homogeneous 
polynomials of degrees n, n and m, respectively, with n < m. We prove that this 
system has at most one limit cycle and that when it exists it can be explicitly 
found. Then we study a particular case, with n = 3 and m = 4. We prove that 
this quintic polynomial system has an explicit limit cycle which is not algebraic. 
To our knowledge, there are no such type of examples in the literature. 
The method that we introduce to prove that this limit cycle is not algebraic can 
be also used to detect algebraic solutions for other families of polynomial vector 
fields or for probing the absence of such type of solutions. 



1. Introduction and Main Results 

Examples of planar polynomial vector fields having explicit algebraic limit cycles 
appear in most textbooks of ordinary differential equations. One of the simplest 
examples is the one of a cubic system that in polar coordinates writes as f = 
r(l — r 2 ), 9 = 1, see for instance [B]. On the other hand, it seems intuitively 
clear that "most" limit cycles of planar polynomial vector fields have to be non- 
algebraic. Nevertheless, until 1995 it was not proved that the limit cycle of the 
van der Pol equation is not algebraic, see [7j. 

The goal of this paper is to give a planar polynomial vector field for which we 
can get an explicit limit cycle which is not algebraic. As far as we know there are 
no examples of this situation in the literature. 

Recall that a real or complex polynomial F(x, y) is an algebraic solution of a 
real polynomial system (x,y) = (X(x,y),Y(x,y)) if 

9F{x,y) dF(x,y) 

— — X(x, y) + — — — Y(x, y) = K(x, y)F(x, y), (1) 

for some polynomial K(x, y), called the cof actor of F. Notice that when F(x, y) is 
real, the curve F(x, y) = is invariant under the flow of the differential equation. 
Observe also that the degree of the cofactor is one less than the degree of the 
vector field. A limit cycle is called algebraic if it is an oval of a real algebraic 
solution. 

Our main result is: 
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Theorem 1.1. The planar differential system 

x = — (x — y)(x 2 — xy + y 2 )+x(2x 4: + 2x 2 y 2 + y 4: j, 
y = -(x + y){2x 2 - xy + 2y 2 ) + y(2x A + 2x 2 y 2 + y 4 ), ^ 

has exactly one limit cycle which is hyperbolic and non- algebraic. In polar coordi- 
nates, this limit cycle is 



4/_\ , i \ -V2 



r = e ^- e ( a + 2 / CQS * + 1 ^)-^ 

cos^(s) + 1 



where i>{6) = 2 f° ^j^ds and 



I _ e 4n-6V27T J Q cOS 2 (s) + 1 

The main steps of the proof are: 

(i) We consider the family of systems, 

x = P n (x, y) + xR m (x, y), 



e W{s)-2s ds w L1990 3. 



(3) 

y = Qn{x,y) + yR m (x,y), 

where P n (x,y), Q n (x,y) and R m (x,y) are homogeneous polynomials of de- 
grees n, n and m, respectively, with n < m. In Section El we prove that all of 
them have at most one limit cycle and that when it exists it can be explicitly 
found. 

(ii) We then consider a concrete system of the form (jSJ) with n = 3 and m = 4, 
having an explicit, unique and hyperbolic limit cycle, see Section HJ 

(iii) Finally, we use the method developed in Section for studying the algebraic 
solutions of the fixed system. Once all of them are found, none of which 
has ovals, we can prove that the limit cycle is non- algebraic, as shown in 
Section El 

In Section 01 we study system (jHJ) with n = 1. In this case we prove that when 
the limit cycle exists it is always algebraic. This is a short preliminary that we 
include for the sake of completeness and that helps to understand how we have 
arrived to the final example studied in our main theorem. 

We end this introduction by showing a system of the form ©, almost equal to 
the one studied in Theorem ll.il having also a unique hyperbolic limit cycle, but 
algebraic. Concretely, it is easy to check that the system 

x = —(x — y)(x 2 — xy + y 2 )+x(x 4: + 3x 2 y 2 + 2y 4: ), 
y = -{x + y){2x 2 - xy + 2y 2 ) +y(x 4 + 3x 2 y 2 + 2y 4 ) ^' 

has the algebraic limit cycle 1— x 2 — y 2 = 0. The results when n — 1 and systems (J2J) 
and (j3J) illustrate that the problem of distinguishing whether the limit cycle for 
system © is algebraic or not is quite hidden in the coefficients of the system. 

2. Systems with explicit limit cycles 
This section is devoted to study system ©. In polar coordinates it writes as 

r = f(9)r n + h(6)r m+1 , 

= g^ey™- 1 
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where 

f{6) = cos0P n (cos0,sin0) + sin 9Q n (cos9, sin 9), 
g{9) = cos9Q n (cos9,sm9) — sm9P n (cos9,sm9), 
h{9) = R m (cos9,sm9). 

Theorem 2.1. System has at most one limit cycle. When it exists it is 
hyperbolic and in polar coordinates it writes as 



cxp 



, Hs) ( f s f(w) , . , 
a + ——- exp — / — - — -aw as 



o g( s ) \ Jo g(w) 



lo 9{ s ) 

where a = ABj (1 — A), being 



g(s) J " Jo g(s) \ Jo g(w) 

Proof. Consider the expression of system (JHl) in polar coordinates, i.e. system (0). 
If g(9) vanishes for some 9 = 9* then it has {9 = 9*} as an invariant straight line. 
From the uniqueness of solutions we get that system Q has no limit cycles. If 
g(9) then we can write the system as 



d9 g(9) g(9) 

which is a Bernoulli equation. By introducing the standard change of variables 
p = r n-m-i we obtain the linear equation 

T9 = {n - m - 1) W + {n - m - 1) WY (6) 

Notice that system (jSJ) has a periodic orbit if and only if equation © has a strictly 
positive 27r periodic solution. 

The general solution of equation (jHJ), with initial condition p(0) = po, is 

p{9- p ) = exp ( f ^\ds)(p + T^expf- f f -^\dw ) ds) := G(9, p ). 

\Jo g(s) A Jo g( s ) V Jo g(w) J J 

(7) 

The condition that the solution starting at p = p is periodic reads as A(p + B) = 
Po- Hence, if A = 1 and B = 0, system (jHJ) has a continuum of periodic orbits, 
otherwise it has at most the solution starting at po = AB/(1 — A) := a. 

In order to prove the hyperbolicity of the limit cycle notice that the Poincare 

return map is II(po) = p(2ir;p ). Thus U'(po) = exp (j^ j^dsj = A ^ 1 for 

all po, and in particular we get that the limit cycle is hyperbolic, whenever it 
exists. □ 

^From the proof of the above theorem we also get the following remark. 

Remark 2.2. When system |3J) has a limit cycle it can be written in the form 
F(r, 9) := r m ~ n+l G(9, a) — 1 = 0. As we will see, the function F(r, 9) can be alge- 
braic or not in cartesian coordinates, depending on the concrete system considered. 
In any case the expression given in (QJ) can be also extended to non-algebraic func- 
tions F, and in this case the cof actor K is not necessarily a polynomial. Curiously 
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enough, independently of the algebraiticity of F its corresponding cof actor satis- 
fying (Op is always the polynomial K(x,y) = (m — n + l)R m (x,y). Examples of 
non-algebraic solutions having a polynomial cof actor have been already given in the 
literature, see for instance |2j, jlj. The example presented in Theorem M . 1\ provides 
a non-algebraic limit cycle having a polynomial cofactor. 

The next remark explains why our proof that the limit cycle of system (J2J is not 
algebraic does not use the explicit expression of the limit cycle. On the contrary, in 
Sectional we develop a method for studying all the algebraic solutions of a system 
having at least a solution of the form y — a(x) = 0, where ct{x) is a rational 
function, and we apply this method to determine all the algebraic solutions of 
system (J2J). 

Remark 2.3. Although for system we know that the expression of the limit 
cycle is r m ~ n+1 G(8, a) — 1 = ; it is not an easy task to elucidate whether this 
curve is algebraic or not in cartesian coordinates. As a example of this difficulty 
we recall the Filiptsov's example, see 0: 

(x = 6(l + a)x + 2y - Q(2 + a)x 2 + 12xy, 

\y = 15(1 + a)y + 3a(l + a)x 2 - 2(9 + 5a)xy + 16y 2 , 

which has the algebraic solution 3(1 + a) (ax 2 + y 2 ) + 2y 2 (2y — 3(1 + a)x) . This al- 
gebraic solution contains a limit cycle for < a < 3/13. For the sake of simplicity 
we fix a = 1/6. For this value the limit cycle is rG(8) — 1 = 0, where 

7(sin 4 g-2sin 2 fl + 1) 

~ 6sin#(-17sin 2 # + 42sin0cos#-7±2 v /^) 

and 

<p(6) = 60 sin 4 9 - 357 sin 3 9 cos 9 + 84 sin 2 9 + 441 sin 2 9 cos 2 9 - 147 sin 9 cos 9. 

Note that is not easy at all to realize that the expression rG(9) — 1 = corresponds 
to a polynomial in cartesian coordinates. 

It is not difficult to see that system (J3J) has always algebraic solutions. 

Lemma 2.4. System (Gj) has F(x,y) = yP n (x,y) — xQ n (x,y) as an algebraic 
solution with cofactor (n + l)R m + div(P n ,Q n ). Notice that it is formed by a 
product of (complex or real) invariant straight lines through the origin. 

Proof. By using the homogeneity of P n and Q n , we know from the Euler's formula 
that nP n (x, y) = xff + y^ and nQ n (x, y) = x^ + y^. Thus 

( dP n dQ n \ ( dP n dQ n \ 

\ V ~dx ~^ n ~ x ~dx~ ) ^ % ' 1 y ~dy~ ~ X ~dy~ ) ^ V ' 

(( xn B dP n dQ n \ 

□ 

Lemma 2.5. Let F(x,y) be an algebraic solution of degree i of the system 

(x = P(x,y) + xR m (x,y), 
\ V = Q(x,y) + yR m (x,y), 
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where P(x,y) and Q(x,y) are polynomials of degree less or equal than n and 
R m (x,y) is a homogeneous polynomial of degree m, with n < m. Thus the homo- 
geneous part of maximum degree of its cof actor is £R m (x, y). 

Proof. Since F is an algebraic solution of system (j3J) we know that 
OF dF 

— (p + X R m ) +^j(Q + yRm) = kf, 

where K is the cofactor of F. 

Denote by F^(x, y) and by K m (x, y) the homogeneous parts of maximum degree 
of F(x,y) and K(x,y), respectively. By using the homogeneity of F e we know, 
from the Euler's formula, that IFi = + V^§f- By equating the higher degree 
terms in the above equation we obtain 

dF £ dF e 
£F £ R m = ——xR m + —-yR m = K m F t . 
ox ay 

Thus K m (x,y) = £R m (x,y) as we wanted to prove. □ 

Next lemma collects some easy remarks on the structure of the cofactors. 
Lemma 2.6. Let 

x = X(x,y), 
V = Y(x,y), 
be a real planar polynomial system. The following holds: 

(i) If it has a complex algebraic solution, then it also has a real algebraic solution. 

(ii) Assume the vector field satisfies 

(X(—x, —y), Y(—x, —y)) = (— y), Y(x, y)), being s either or 1. 

Then if it has a real algebraic solution then it has another real algebraic 
solution with cofactor K satisfying K(—x, —y) = (—l) s+1 K(x,y). 

Proof. The proof is elementary. 

□ 

Finally we give an integrating factor for system (JHJ). 

Lemma 2.7. Consider system and define 

V(x,y) = (r m - n+1 G(6,a) - l)(yP n (x, y) - xQ n (x, y)), 

where G(9,po) is the function given in ffl) and po = a is the value for which this 
function is 2n— periodic. Then, whenever it is defined, 1/V(x,y) is an integrating 
factor of the system and we call V(x, y) an inverse integrating factor. 

Proof. We use the following formula: let Fx and F 2 be two solutions of (x, y) = 
(X(x,y),Y(x,y)) with cofactors K\ and K 2 , respectively. Thus 

div (^r^) = (div(X, Y) - (K, + K 2 )) . 

We remark that the above formula, taking a denominator of the form [Ji^, for 
some real or complex constants cti, is indeed the key point of the Darboux theory 
of integrability, see 0. Take F\(x,y) = yP n (x,y) — xQ n (x,y) and F 2 (x,y) = 
r m ~ n+l G(8, a) — 1. By using Lemma 12.41 and Remark 12.21 we know that their 
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associated cofactors areit being non-Liouvillian Ki(x,y) — (n + l)R m (x,y) + 
dxv(P n (x,y),Q n (x,y)) and K 2 (x,y) = (m - n + l)R m (x, y), respectively. On 
the other hand, taking the vector field associated to system Q we get 

f)R fin 

div(X, Y) = div(P„, Q n ) + 2R m + x—^ + y—^ = div(P n , Q n ) + (2 + m)R m , 

ox ay 

where we have used again Euler's formula. Collecting all the above results we get 
div ((X, Y) I [FxF-ij] = as we wanted to prove. □ 

Remark 2.8. (i) When we apply the above Lemma to systems (HJ) and O we get 
both non-algebraic and algebraic inverse integrating factors. 

(ii) In [5J it is proved that when 1/V(x,y) is an integrating factor of (x,y) = 
(X(x,y),Y(x,y)) and V(x,y) is defined in the whole plane, all the limit cycles 
of the system are included in the curve V(x, y) = 0. This is the case of sys- 
tem the limit cycle, whenever it exists, it is given by the expression F 2 (x, y) = 
r m - n+1 G(9, a) - 1 = 0. 

(Hi) The equality div ((X, Y) / (F 1 F 2 )) = also holds when instead of F 2 (x, y) = 
r m ~ n+1 G(9, p )\ po=a — 1 we take a different value of p , but in this case F 2 is indeed 
a multivaluated function and the result of jH] can not be applied. 

3. A FAMILY OF SYSTEMS WITH EXPLICIT ALGEBRAIC LIMIT CYCLES 

The existence of limit cycles for a subfamily of system (JHJ) has been studied in 

0. Here we prove that the limit cycle found there is algebraic. 

Proposition 3.1. Consider the system 

x = -y + x (a + R m (x,y)), , g . 
y = x + y(a + R m (x,y)), 

where a is a real parameter and R m (x, y) is a homogeneous polynomial of degree m. 
Then it has only two algebraic invariant curves x 2 + y 2 and H(x, y) = G m (x, y) — 

1, where G{6) = G m (cos9,sm6) satisfies G' + maG + mR m (cos 9, sin 6) = 0. 
Furthermore, when the limit cycle exists, m is even and H(x, y) contains a real 
oval which is the limit cycle of the system. 

Proof. iFrom Lemma 12.41 we can see that x 2 + y 2 is an algebraic solution with co- 
factor K(x, y) = 2a + 2R m (x, y). Now, we study other possible algebraic solutions. 
Write the Fourier expansion of R m : 

k=m 

R m (cos9, sin#) = CkC ke \ where C/T = c_ fe 6C,cj = when k ^ m(mod 2). 

k=—m 

Following the steps of the proof of Theorem 12.11 we obtain that, in polar coordi- 
nates, the solution of (JHJ) starting at r = ro when 9 = can be written as 

k=m 



r 



-m | —m 





V — — e~ mae + G m (cos 9, sin i 

' J tci -1- rnn I 



ki + ma 



or 



k=m 

+ m V — — ) r m e- ma9 + G m (rcos^,rsini 

^— ' ki + ma I 
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where G m (x, y) is the homogeneous polynomial of degree m defined by its Fourier 
expansion as 

k=m 

G m (cos9,sm9) := -m V e k9i , 

' fcz + ma 

fc=— m 

and G(6) = G m (cos9, sm9) satisfies G' + maG + mR m (cos9, sm9) = 0. 

By using the above expression we get that the only algebraic solution of sys- 
tem (EI) is the one that satisfies 



k—m 



r-™ + m y — = 0. 

^— i hi -\- m n 



ki + ma 

K=—m 

Moreover, it is easy to check that the cofactor of this algebraic solution, H(x, y) = 
G m (x, y) — 1, is K(x, y) = mR m (x, y), see also the proof of Lemma l2~71 

Notice that a necessary and sufficient condition for the existence of a real alge- 
braic solution is that 

k=m 

E 7-^<0. 
^— ' fez + ma 

k=—m 

Finally, the limit cycle exists when G m (cos 9, sin 9) > for 9 e [0, 27r]. This only 
can happen when m is even, see also [3]. □ 

4. The examples 

In this section we will prove that the system given in Theorem ll.ll has an explicit 
limit cycle. That the limit cycle is not algebraic will be proved in Sectional 
The system (J2J) is 

x = — ( y x — y)(x 2 —xy + y 2 )+x(2x 4 + 2x 2 y 2 + y 4 ), 
y = -{x + y)[2x 2 - xy + 2y 2 ) + y[2x 4 + 2x 2 y 2 + y 4 ), 



and in polar coordinates it can be written as 



r = (cos 4 9 + l)r 5 + (cos 2 #-2)r 3 , 
= -(cos 2 9 + l)r 2 . 

By following the same steps than in the proof of Theorem 12.11 we introduce the 
change of variables r = 1/ v /p, obtaining: 

,_ cos 2 #-2 cos 4 # + l 
9 ~ cos 2 # + l P + cos 2 9 + 1 
The solution satisfying that p = p > when 9 = is 

p{0- p ) = e~^ 2d (p + 2 f C °*l\ S \ + \ ^- 2 °ds) > 0, 

where ^) = 2jJ^>. 

The initial condition of the limit cycle is given by the equation p(2tt) = p(0) = 
Pq. Hence 

ro 1 - e^-e^ J cos 2 (s) + 1 



<s 
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This value can be can be computed numerically, giving p* m 1.1990. The intersec- 
tion of the limit cycle with the OX + axis is the point having rg = l/^Po ~ 0.9132. 

Since the Poincare return map is n(p ) = p(27r; p ) we have n'(p ) = e 1 - 4-6 ^ 77 < 1 
for all p an d < 0, we get that the limit cycle of system (J2J) is hyperbolic and 
unstable. 

5. A METHOD FOR STUDYING THE EXISTENCE OF ALGEBRAIC SOLUTIONS 

Let F(x,y), K(x,y), X(x,y) and Y(x,y) be real analytic functions such that 
OF OF 

— X(x, y) + -q-Y(x, y) = K{x, y)F{x, y). (9) 

Thus it is clear that the set {(x, y) G 1R 2 : F(x,y) = 0} is formed by solutions of 
the system 

x = X(x,y), 



V = Y(x,y). ^ 

Fixed an analytic solution of fllOl) of the form y = a(x), we can consider the 
following Taylor expansions in z, 

F(x,z + a(x)) = F Q (x) + zFt(x) + z 2 F 2 {x) + . . . , 

K(x,z + a(x)) = K Q (x) + zKx{x) + z 2 K 2 (x) + . . . , 

X(x,z + a(x)) = X Q (x) + zXt(x) + z 2 X 2 (x) + . . . , 

Y(x,z + a(x)) = Y (x) +zYi(x) + z 2 Y 2 {x) + ... 



Notice that a'(x) = Y (x)/X (x). Then equation Q can be written as 

oo / k \ 

\ Y.( Xk -^ x ) F '^ + i iY ^+i - ia'(x)X k ^ i+1 - K k -i)Fi(xj) )z k = 0. 

k=0 \i=0 J 

The functions F k (x) can be obtained recurrently from the above relation by solving 
the linear differential equations in Fk(x), obtained vanishing each coefficient in z k . 
In particular, for k = and k — 1, we get 

X (x)F^x)-K (x)F Q (x) = 0, 
X (x)F[(x) + (Y 1 (x)-a'(x)X 1 (x)-K (x))F 1 (x) + X 1 (x)F^x)-K 1 (x)F (x)=0. 
We obtain F (x) = Co exp ( J Q X ^44 ds) , where Co is an arbitrary constant, and 
similarly we could get Fi(x). 

When a(x) is a rational function and F(x,y), K(x,y), X(x,y) and Y(x,y) are 
polynomials, with real or complex coefficients, the linear differential equations 
for each F k (x) described in the above algorithm give us a collection of necessary 
conditions for the existence of an algebraic solution F(x,y). The conditions are 
that, for each k, the functions Fk(x) must be polynomials. For instance, for k = 0, 
the first necessary condition is that the primitive of the rational function 

K (x) K(x,a(x)) 
Xq(x) X(x,a(x)) 

must be a linear combination of logarithms of polynomials. Furthermore the 
coefficients of the logarithms have to be natural numbers. 
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The necessary conditions obtained for the existence of algebraic solutions restrict 
the possible cofactors of F. These restrictions give the key for searching the possible 
algebraic solutions of system (fTUj). see Remark 15.21 and Sectional As we will see, 
in our case we only need to apply the described method for k = but we remark 
that in other situations, by using it for bigger k, it can give more information 
about the existence or non existence of algebraic solutions. 

Remark 5.1. Notice that the above method can only be applied when the candidate 
F(x, y) to be an algebraic solution of system M0\) does not contain the factor 
V - oc(x). 

Remark 5.2. Assume that system il(J\) is fixed and it is polynomial. Notice that 
the equation (0) that gives the possible set of algebraic solutions of system MU\) is 
equivalent to a set of quadratic equations where the unknowns are the coefficients 
of F and the coefficients of K. In general, it is very hard to solve this system of 
equations, even by using algebraic manipulators. On the other hand, the method 
developed in this section imposes restrictions on the cof actor K for the existence 
of F. Ideally, if K is totally known, the system to be solved is be linear and so the 
problem of knowing the existence or not of algebraic solutions of a given degree 
would be a much easier task. In any case, any information on K makes the problem 
simpler. 



6. A NON-ALGEBRAIC LIMIT CYCLE 

This section is devoted to prove that the limit cycle of system (J2J) is not algebraic. 
Indeed, by using the method introduced in Section El we will prove that the only 
algebraic solutions of the system are the ones given in Lemma f2. 41 The algebraic 
solution given by this lemma is yP n — xQ n = {2x 2 + y 2 ){x 2 + y 2 ). Concretely, the 
curves x 2 + y 2 = and 2x 2 + y 2 = have the cofactors 2(2x 4 + 2x 2 y 2 + y 4 — x 2 — 2y 2 ) 
and 2(2x 4 + 2x 2 y 2 + y 4 — x 2 + xy — 2y 2 \ respectively. These two curves coincide 
with the four complex lines y = ±ix and y = ±\/2ix. 

As we will see, the first step (k = 0) of the method developed in Section 
applied to each one of the four complex lines, y = ±ix and y = ±\/2ix, will give 
enough restrictions to prove that the only algebraic solutions of system (j2J) are the 
ones described above. 

Assume that the differential system has a real or complex algebraic solution F 
and that it does not contain none of the given four lines as a factor. By using 
Lemma \2 .61 it is not restrictive to assume that F is real and that its cofactor is an 
even function, i.e. K(—x, —y) = K(x,y). 

Since the degree of the vector field © is 5 we know that the degree of K(x, y) 
is at most 4. By the above restrictions on K(x, y) and by using also Lemma 12.51 
we can write it as the real polynomial 

K(x, y) = a o + a 2 ox 2 + a u xy + a 02 y 2 + i{2x 4 + 2x 2 y 2 + y 4 ), 

where I is the degree of the corresponding algebraic curve F(x, y) = 0. 
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We apply the first step of our method, i.e. we take k = 0. By considering the 
cases a(x) = ±ix we obtain 

K(x,a(x)\ dx = f K(x,±ix) = f Kg) = 



X(x,a(x)) J X(x,±ix) J X (x 

a o(— 1 ± i) 1 , , ., w -\ f \ 
^ V ^(^o + an - a 02 ± H~ a 2o + a u + a 02 + a 00 )) log(ar) + 

^(-a 20 - On + a 02 + It ± i(a 20 - a u - a 02 - a 00 )) log(2 + 2x 2 + x 4 ) + 

o 

^(a 2 o - au - «02 - aoo ± i(a 2 o + an - a 02 - 2^)) arctan(x 2 + 1). 

By forcing F(x, a(x)) = F(x, ±ix) = F (x) = C exp (J K (x, ±ix)/X (x, ±ix) dx) 
to be a polynomial, with Co an arbitrary constant, we obtain a first set of necessary 
conditions: 

a 2 o — On — a-02 — «oo = 0, 
«2o + on — o,02 — 2£ =0, 
a 00 = 0. 

The same computations can be done for the other pair of algebraic solutions, 
y = ±\^2ix, that is 

*(»,<«(»)) (fc= r^,±^) ix= t^ dx= 



X(x,a(x)) " J X(a;, ±zV2a;) ' J X (x 
/ "° + ^( 3a 20 - 6a 02 - 2a 00 ± 3V2ia u ) log(x) 



6x 2 9 

— (-3a 20 + 6a 02 + 9£ + 2a 00 =p 3V2ia u ) log(3 + 2a; 2 ). 

lo 

As in the previous case, we obtain a second set of necessary conditions: 

an = 0, 
a 00 = 0. 

Collecting all the obtained equations we get that the degree of the invariant alge- 
braic curve is i = 0, or in other words that such a curve does not exist. 

Proof of Theorem 11.11 The proof of the theorem follows from the results of 
Sections H Hand El 

Remark 6.1. $From (HI Thm 1], if a planar system has an explicit non-algebraic 
solution which is in the zero level set of a Liouvillian function then it has a Dar- 
boux integrating factor and therefore the whole system is integrable by quadratures. 
Notice that this is the case for system it has a non- algebraic Liouvillian limit 
cycle and it can be transformed into a Bernoulli equation. Consequently, if we 
would like to have an explicit non-algebraic limit cycle for a planar system, which 
is not integrable by quadratures, we should look for a limit cycle given by a non- 
Liouvillian function. 
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